Cosmic Ray Studies with the Fermi Gamma-ray Space Telescope 

Large Area Telescope 


D. J. Thompson*, L Baldini b , Y. Uchiyama 0 

0 NASA Goddard Space Flight Center, Greenbelt, MD 20771 USA 
b Istituto Nazionale di Fisica Nucleate, Sezione di Pisa, 1-56127 Pisa, Italy 
C SLAC National Accelerator Laboratory, 2575 Sand Hill Road M/S 29, Menlo Park, CA 94025 USA 


Abstract 

The Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope provides both direct 
and indirect measurements of Galactic cosmic rays (CR). The LAT high-statistics observations 
of the 7 GeV - 1 TeV electron plus positron spectrum and limits on spatial anisotropy constrain 
models for this cosmic-ray component. On a Galactic scale, the LAT observations indicate that 
cosmic-ray sources may be more plentiful in the outer Galaxy than expected or that the scale 
height of the cosmic-ray diffusive halo is larger than conventional models. Production of cosmic 
rays in supernova remnants (SNR) is supported by the LAT y-ray studies of several of these, both 
young SNR and those interacting with molecular clouds. 
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1. Introduction 

Among the earliest stimuli for studying cosmic y rays was the recognition by Hayakawa [1] 
and Morrison [2] that cosmic ray particles interacting with interstellar matter would inevitably 
produce high-energy y rays. The first high-energy y-ray sky map by OSO-3 [3] confirmed this 
expectation. Gamma-ray and cosmic-ray studies, both observational and theoretical, have been 
closely linked ever since, with principal results for photon energies above 30 MeV coming from 
satellite detectors (SAS-2 [4], COS-B [5] , EGRET/Compton Gamma Ray Observatory [6]) and 
ground-based TeV telescopes (for a recent review, see [7]). 

While sharing many experimental techniques, cosmic-ray and y-ray studies offer comple- 
mentary approaches. Cosmic-ray detectors measure the particles directly in the solar neighbor- 
hood, including detailed spectral and composition analysis. Observations with y-ray telescopes 
study cosmic-ray interactions at a distance, providing indirect insight about the particles. The 
y-ray observations have the potential to trace sources of cosmic rays, since the photons are not 
deflected by Galactic magnetic fields as the charged particles are. 

The present review concentrates on results from the Large Area Telescope on the Fermi 
Gamma-ray Space Telescope satellite, which was launched on 11 June 2008. The Fermi LAT 
offers three lines of research related to cosmic rays: 

1. Because the LAT is a particle detector, it can make direct measurements of cosmic rays, 
especially cosmic-ray electrons (CRE). 

2. The large-scale y-ray properties of our own and other galaxies offer some insight into the 
distribution of cosmic rays and their sources. 
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3. The LAT measures the y-ray properties of likely acceleration/interaction sites of cosmic 
rays such as supernova remnants. 


2. The Fermi Large Area Telescope 

The LAT is a y- ray telescope [8] , designed to distinguish photons in the energy range 20 MeV 
to more than 300 GeV from the high background of energetic charged particles at the 565 km 
altitude orbit of the Fermi satellite. The tracking section has 36 layers of silicon strip detectors to 
record the tracks of charged particles, interleaved with 16 layers of tungsten foil (12 thin layers, 
0.03 radiation length, at the top or Front of the instrument, followed by 4 thick layers, 0.18 
radiation length, in the Back section) to promote y-ray conversion into electron-positron pairs. 
Beneath the tracker is a calorimeter comprised of an 8-layer array of Csl crystals (1.08 radiation 
length per layer) to determine the y-ray energy. The tracker is surrounded by segmented charged- 
particle anticoincidence detectors (ACD; plastic scintillators with photomultiplier tubes) to reject 
cosmic-ray background events. The LAT’s high sensitivity stems from a large peak effective 
area (~8000 cm 2 ), wide field of view (~2.4 sr), good background rejection, superior angular 
resolution (68% containment angle ~0.6° at 1 GeV for the Front section and about a factor of 2 
larger for the Back section), and excellent observing efficiency (keeping the sky in the field of 
view with scanning observations). In scanning mode, the LAT observes the entire sky every 2 
orbits (about 3 hours). Pre-launch predictions of the instrument performance are described in [8] . 

In addition to being a y-ray telescope, the LAT is intrinsically an electron detector. One 
of the biggest analysis challenges for electron studies is the separation of CRE from the much 
more abundant hadrons (primarily protons). Key to the background rejection is the LAT capabil- 
ity to discriminate between electromagnetic and hadronic showers, based on the different event 
topology in the three subsystems. This is largely in common with the photon analysis, the most 
prominent difference being the use of the ACD. 

Unlike previous generations of pair conversion telescopes, the LAT basically triggers on all 
the charged particles crossing the active volume, with no harsh photon selection built into the 
hardware trigger. This drastic change of approach is largely due to the use of silicon detectors, 
allowing precise tracking with essentially no detector-induced dead time. Some event filtering is 
necessary prior to downlink in order to fit the data volume into the allocated telemetry bandwidth. 
The onboard filter is really tailored to preserve the y-ray efficiency by discarding events that are 
unlikely to be y rays. The rate of events depositing more than 20 GeV into the calorimeter, 
however, is small enough that we can afford effectively to disengage the filter above that energy, 
transmitting to the ground high-energy charged particle events along with the y rays. 

3. Direct Cosmic Ray Measurements 

A century after the pioneering cosmic-ray studies by Domenico Pacini and Victor Francis 
Hess we have a huge body of knowledge about these energetic particles. In the LAT energy 
range, many instruments have measured the energy spectra of CR particles and their compo- 
sition, including not only protons and electrons, but also heavier nuclei. They can be roughly 
divided into inclusive measurements (i.e. p and e~ + e~ spectra) and measurements of the charge 
composition (i.e. p/(p + p) and e + /(e~ + e + )), providing complementary information on the ori- 
gin, acceleration and propagation of cosmic rays. This dichotomy has a parallel in the variety of 
experimental techniques routinely exploited since the mid-1960’s, in that cosmic-ray instruments 
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generally fall in one of the two categories: calorimetric experiments (that cannot distinguish, per 
se, the charge sign) and magnetic spectrometers. 
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Figure 1: Compilation of the characteristics (in terras of geometric factor and collection time) of some balloon and 
satellite cosmic-ray and y-ray detectors. Care has been taken in extracting the relevant information from the original 
publications in a consistent fashion. We believe that the information in the plot is accurate within a factor of ® 2, 
which is adequate for the illustrative purpose of the plot itself. Filled markers indicate space-based experiments, while 
open markers indicate balloon experiments; conversely squares indicate calorimetric experiments, while circles indicate 
magnetic spectrometers. 


Figure 1 is a partial compilation of such instruments, classified by acceptance 1 and effective 
observation time. It is not a trivial task to give an exhaustive and unbiased account of the huge 
experimental effort that turned the first observations into precision science. In many cases the 
geometry factor is energy dependent, and some modem instruments (particularly the magnetic 
spectrometers) can operate in different configurations with different acceptances. In addition, 
each detector is sensitive in a different energy band (for any given observable), which makes a fair 
comparison extremely difficult. Despite these caveats, Figure 1 is accurate enough to illustrate 
some interesting aspects of the historical development of cosmic-ray measurements: calorimetric 
experiments feature, in general, a larger acceptance than magnetic spectrometers; and space 
experiments have the obvious advantage of a relatively longer collection time. Figure 1 also puts 
Fermi into the context of this development: the LAT’s unique combination of large acceptance 
and long observation time results in an unrivaled exposure factor 2 exceeding 10 8 m z sr s, as 
indicated by the diagonal gray lines. 


‘The acceptance, or geometry factor (Gy), is the integral of the effective area over the instrument field of view and 
constitutes the proportionality factor between the event rate in the detector and the incident flux. 

2 Tbe exposure factor is the product of the geometry factor and the effective observation time and is directly related to 
the total number of events recorded by a given instrument. 
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4. Direct CR measurements with Fermi 


High-energy cosmic-ray electrons and positrons constitute a peculiar component of the cos- 
mic radiation in that, unlike protons and heavier nuclei, they rapidly lose energy by synchrotron 
radiation on Galactic magnetic fields and by inverse Compton (IC) scattering on the interstellar 
radiation field. In the energy range of interest for the LAT— i.e. above a few GeV, where losses 
dominate over escape— the average lifetime for electrons against synchrotron and IC losses is 
Tf.{E) oc E ~ l , whereas the corresponding escape time for protons scales with the inverse of the 
diffusion coefficient: t p (E) oc E~ s (with <5 w 0.3-0.6). This implies that the all-electron spec- 
trum at Earth is significantly steeper than the proton spectrum. More importantly, at energies 
above a few hundred GeV the majority of the electron and positron flux must be contributed by 
sources closer than a few hundred pc. In other words, high-energy CRE really probe cosmic-ray 
production and propagation in nearby or local Galactic space. 

The other piece of conventional wisdom is that unlike electrons, which are accelerated in 
CR sources, the vast majority of positrons arise as secondary products of cosmic-ray interactions 
in our Galaxy. In the absence of primary sources of positrons it is straightforward to demon- 
strate that in the simplest scenarios the positron fraction is expected to decrease with energy like 
E" e , e.g. [9]. As we shall see in a moment, this expectation is at odds with the most recent 
measurements [10] . 

Fermi began its operations at about the same time two important measurements [11,10] were 
published that challenge the simplest models of cosmic-ray production and propagation. These 
reports stirred an enormous excitement in the community, mainly in connection with the possi- 
bility of detecting signs of new physics in the context of indirect dark matter searches. In 2009 
the Fermi LAT collaboration published the first systematics-limited measurement of the inclusive 
spectrum of cosmic-ray electrons and positrons [12] between 20 GeV and 1 TeV. Though harder 
than previously thought, the Fermi spectrum does not show evidence of any prominent feature, 
such as the one reported by the ATTC experiment [11]. This conclusion was later independently 
confirmed by the H.E.S.S. [13] and Pamela [14] experiments, as well as by a dedicated analysis 
of the Fermi data aimed at the highest end of the energy range [15]. 

The Fermi measurement— especially when the low-energy extension of the spectrum pub- 
lished in [15] is taken into account— is in tension with a single power-law model, though this 
cannot be definitely ruled out within the systematic uncertainties. This tension has motivated a 
number of papers in which a consistent interpretation of the Fermi inclusive spectrum and the 
rise in the positron fraction above * 10 GeV measured by Pamela [10] is attempted through the 
introduction of an additional nearby source of high-energy electrons and positrons (for a more 
detailed discussion see [17] and references therein). It is fair to say that, to date, the issue is not 
settled, and there is no firm consensus on whether observations really indicate the existence or 
nature of this extra component. Indeed there are alternative scenarios (see for instance [18, 19]) 
that might explain the observations without invoking any new (astrophysical or exotic) source of 
high-energy electrons and positrons. 

For the future, with more than 5000 electrons above 1 TeV per year crossing the LAT, Fermi 
can hope to provide an independent confirmation of the spectral cutoff measured by H.E.S .S. [20] 
in this energy range. Even with a selection efficiency on the order of 10% (in line with the 
figures of merit of the current analysis), the statistics will not be an issue. At the same time, the 
collaboration is actively working on the energy measurement, which is the most challenging part 
of the analysis. 

Without an onboard magnet, Fermi LAT cannot directly distinguish the charge sign. How- 
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Figure 2: Cosmic-ray electron plus positron spectrum measured by Fermi in one year of observation, adapted from [IS]. 
Some of the most recent measurements are shown for comparison. The dashed line is a pr e-Fermi model, shown 
for reference. Many of the data points are taken from the cosmic-ray database maintained by A. W. Strong and 
I. V. Moskalenko [16], 
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ever the Earth’s magnetic field creates natural shadow regions from which particular charges 
are forbidden— simply because their paths are blocked by the Earth— and can be effectively ex- 
ploited in order to separate positrons from electrons. The idea of using the geomagnetic field to 
study the charge asymmetry in the leptonic component of cosmic rays is not new. In fact the pos- 
sibility of using the east-west asymmetry in the geomagnetic cutoff rigidity has been explored in 
several balloon experiments (see [21] and references therein for more information). In addition 
to its large field of view and its excellent directional capabilities, Fermi can take advantage of the 
accurate models of the Earth’s magnetic field [22] and the particle tracing computer programs 
that are nowadays readily available. One of the main challenges in the analysis is the fact that, as 
the energy increases, the shadow regions become smaller and smaller. Coupled with the falling 
spectrum, this effect in practice limits the highest accessible energy to a 200 GeV. In addition 
to that, these regions are located near the Earth’s horizon, which is outside or at the limits of the 
LAT field of view in the nominal sky-survey operations. Preliminary results of this analysis are 
in agreement with the positron fraction measured by Pamela and will be the subject of a forth- 
coming paper. With the AMS detector now installed on the International Space Station, new data 
in the unexplored high-energy range will be available in the near future. 

The electron component of cosmic rays is not fully described by the inclusive energy spec- 
trum and the positron fraction: the degree of anisotropy in the arrival directions is an additional 
key observable, especially in a scenario in which the inclusive electron spectrum is dominated 
by one (or a few) nearby sources. Some degree of anisotropy is indeed expected even in a plain 
diffusive scenario in which the sources are continuously distributed, just due to the fact that the 
Earth is not at the center of the diffusion halo. It is natural to expect that the flux anisotropy can 
be much larger if a single (or a few) nearby sources dominate the CRE spectrum at high energy. 



Figure 3: Fermi LAT 95% CL upper limits for the integral dipole anisotropy 6 as a function of the minimum energy 
(adapted from [23]). For comparison the expected anisotropy levels for two pulsars (Vela and Monogem) that are plau- 
sible candidate sources of high-energy electrons are shown. The values of the model parameters are described in [23); 
though subject to large uncertainties, the values are physically motivated and respect the relevant constraints (the first of 
which is the spectrum measured by Fermi). The (much lower) expected anisotropy for a continuous source distribution 
throughout the Galaxy (calculated with GAI.PROP [24]) is also shown for reference. 

Key to these measurements are Fermi’s large geometric factor, its excellent directional capa- 
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bilities and its characteristic of an all-sky monitor providing uniform exposure 3 . 

Fermi has published the most stringent upper limits to the anisotropy of the CRE flux, based 
on the first year of data [23], The paper discusses all aspects of the analysis: the pixelization 
of the sky, the pixel-to-pixel fluctuations, the search for anisotropies on different angular scales 
(from 10° to 90°) and a full spherical harmonic analysis. In this review we concentrate on a single 
aspect, namely the upper limits on a possible dipole anisotropy, which is physically relevant for 
a single source dominating the spectrum (in which case the source defines one of the poles of the 
dipole). 

Some of the basic results are shown in Figure 3. Though entering a detailed discussion of 
specific models and related parameter space is beyond the scope of this review, we emphasize 
that Fermi is already starting to constrain physically motivated scenarios and might hope to detect 
a positive signal during the lifetime of the mission (5 to 10 years). If indeed there is a nearby 
source of high-energy electron/positron pairs, this signal would be a powerful tool to shed light 
on its nature. For most of the dark matter scenarios one would expect an excess toward the 
Galactic center (potentially at the same level as the ones shown in the figure) while Monogem 
and Geminga, which are two good pulsar candidates, are roughly in the opposite direction. 

5. Diffuse Galactic Gamma Rays as Probes of Galactic Cosmic Rays 

Unlike the sky at visible wavelengths, the y-ray sky is strongly dominated by diffuse radia- 
tion originating in our Milky Way Galaxy, making the Galactic plane the most striking feature 
in any large-scale y-ray skymap. Although the diffuse Galactic y radiation is brightest along 
the plane, it is seen at all Galactic latitudes. This radiation is largely produced by cosmic-ray 
interactions with the interstellar gas and photon fields through the processes of inelastic nucleon 
scattering, bremsstrahlung, and inverse Compton scattering. Figure 4 shows the Spectral Energy 
Distribution (SED) of this radiation in part of the sky, along with modeled components [25]. 
The largest contributor to this radiation is nucleon-nucleon collisions (basically cosmic ray pro- 
tons and heavier nuclei hitting interstellar hydrogen nuclei) with subsequent decay of rr 0 mesons 
into y rays in the Fermi LAT energy range. The LAT measurements do not confirm the unex- 
pected excess in the few GeV energy range measured by EGRET on the Compton Gamma Ray 
Observatory (CGRO) compared to models. 

Mapping the Galactic y rays offers a way to derive information about the spatial distribution 
of cosmic-ray particles. The modeling necessary to extract such results is not simple because it 
requires knowledge of all the matter and photon distributions with which the cosmic rays interact. 
Models for distribution of cosmic-ray sources and propagation include GALPROP [26] 4 and a 
derivative, DRAGON [27]. Figure 5 shows one result that has been derived from the LAT results: 
an indication that conventional models for CR source distributions underpredict the observations 
in the outer Galaxy. As noted in [28] , the observed y-ray emissivity gradient is inconsistent with 
all available tracers of massive star formation. In this case the LAT results suggest either a large 
scale height for the CR diffusive halo or a flatter distribution of the sources beyond the Solar 
circle, although other possibilities cannot be ruled out, including a large mass of unseen gas in 
the outer Galaxy [29], a non-uniform diffusion coefficient [27], or advective CR transport by a 
Galactic wind [30] . 


3 The exposure is not actually uniform over the sky; there is a non-uniformity at the level of 15%, mainly due to the 
fact that the instrument is not taking data while inside the South Atlantic Anomaly. This is described in [23). 

4 see also http://galpropjlanford.edu/ 
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Figure 4: Spectral Energy Distribution of diffuse y radiation at intermediate Galactic latitudes. Overlaid jt SEDs for 
the component processes of the diffuse emission. The UIB component is unidentified isotropic background, including 
instrumental background. The primary component, shown by the heavier line above the gray band, is decay of n° 
mesons [25], 


Some galaxies not exhibiting apparent nuclear jet activities are also seen as y-ray sources 
that offer information about cosmic rays beyond our own Galaxy. Examples are normal galax- 
ies, including the Large Magellanic Cloud [35], the Small Magellanic Cloud [34], and starburst 
galaxies, such as M82 and NGC 253 [36], GeV y rays in these galaxies come primarily from the 
interactions of cosmic ray hadrons and electrons with interstellar matter and photon fields. Fig- 
ure 6 shows the measured spectrum from the Small Magellanic Cloud. The spectrum is modeled 
predominantly by 7t° decay, implying hadronic processes similar to those in our own Galaxy. The 
detailed analysis of the LMC shows that the CR flux there is less than a third that of the Milky 
Way [35], and the derived CR flux in the SMC is even lower [34] . The y-ray luminosities of these 
normal and starburst galaxies show an approximately linear relationship with the product of the 
supernova rate and the total mass of gas in the galaxies [36] . The suggestion is that the y-ray 
production requires both target material (gas) and CR production (traced by the supernova rate). 
Although this analysis is highly simplified, the trend supports the expectation that the principal 
accelerators of cosmic rays are in some way related to massive-star- formation activity. 

6. Supernova Remnants and Other Candidate Cosmic-Ray Sources 

6.1. SNRs as Cosmic-Ray Sources 

Supernova remnants (SNRs) have long been thought to be the main sources of Galactic 
CRs [37], accelerating CR particles at their expanding shock waves through the diffusive shock 
acceleration (DSA) process [38]. In this standard paradigm, SNRs are responsible for the pro- 
duction of cosmic-ray protons and nuclei up to energies of ~ 3 x 10 15 eV, transferring ~ 10% of 
the explosion kinetic energy into the form of cosmic-ray energy [39]. The observations of y rays 
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Figure 5: CR source distribution from one GALPROP model (solid line), compared with the SNR distribution obtained 
by the £-D relation [31] (see [32] for a critique of this method, however) and that traced by the pulsar distribution L33] 
shown by dotted lines. The thin solid line represents an example of the modified distributions introduced to reproduce 
the emissivity gradient observed by the LAT [28] . 


from SNRs have been viewed as the most promising method of addressing the SNR paradigm for 
the origin of cosmic rays, through the measurements of the 7 r°-decay (or hadronic ) y rays [40], 
Despite the substantia] progress of TeV 7 -ray observations in recent years [7], finding a con- 
vincing case of 7 T° -decay 7 rays remained tantalizingly difficult, mainly because electromagnetic 
radiation processes involving relativistic electrons (the so-called leptonic components) could ex- 
plain the 7 -ray emission as well. It has widely been expected that the Fermi LAT could identify 
the hadronic component of the 7 -ray emission in SNRs. Indeed the observations with the Fermi 
LAT have started to provide an effective means to disentangle the leptonic (IC scattering and 
relativistic bremsstrahlung) and hadronic (fr°-decay) 7 -ray components. 

Gamma-ray observations with the Fermi LAT have so far resulted in the detections of GeV 7 - 
ray emission from more than ten SNRs [41 , 42, 43, 44, 45, 46, 47, 48] . Figure 7 summarizes the 
observed GeV luminosity in 0.1-100 GeV as a function of the diameter squared for the known 
LAT-dctected SNRs. Note that the squared diameter would be regarded as a reasonable indicator 
of SNR ages, as the age of a remnant in the Sedov-Taylor phase is related to its radius as t oc R 2 - 5 . 
The remnants interacting with molecular clouds are indicated by filled circles. They constitute 
the dominant class of the LAT-detected SNRs, and they are generally more luminous in the 
GeV band than young SNRs [49]. The observed 7 -ray luminosity of the SNRs interacting with 
molecular clouds is so high that one of the main leptonic radiation channels, namely IC scattering 
off the Cosmic Microwave Background (CMB) and interstellar radiation fields , becomes uni ikely . 
For an efficiency of ~ 10% for CR proton production and an electron-to-proton ratio of 0.01 , the 
GeV luminosity of the IC emission is calculated to be Lie ~ 10 34 erg s _i using the interstellar 
radiation fields of the solar neighborhood. The observed luminosity well exceeds this value as 
shown in Figure 7. In what follows, we review the results obtained with the LAT in the context 
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Figure 6: SED of y radiation from the Small Magellanic Cloud. Overlaid are SEDs for the component processes of the 
diffuse emission. The primary component is decay of jt° mesons [3d] . 


of the SNR paradigm for CR origin. 

6 2 . Young Supernova Remnants 

The remnants of historical supemovae are generally well studied across the electromagnetic 
spectrum, making them excellent laboratories for studying high-energy phenomena associated 
with supernova shocks, and for testing our current understanding of the DSA theory. The follow- 
ing issues can be addressed by the Fermi LAT observations of young SNRs: 

(i) Acceleration efficiency— About 10% of the mechanical energy of supernova explosions 
must be transferred to the kinetic energy of CR protons in order to maintain the density of Galac- 
tic CRs. Since the detailed mechanisms of how thermal particles are injected into the acceleration 
process and therefore the efficiency of particle injection remain uncertain, the current theory does 
not firmly predict the amount of CR protons in supernova shocks. The amount of CR electrons 
produced at shocks is even more uncertain theoretically; 

(ii) Magnetic field amplification— The DSA process at supernova shocks depends strongly 
on the generation of magnetic field turbulence by CR themselves ahead of the shock wave. It 
has been predicted that the turbulent magnetic field in the strong shocks of young SNRs can be 
significantly amplified by CR current-driven instability [50] and other instabilities. Evidence for 
the magnetic field amplification was provided by high angular resolution observations with the 
Chandra X-ray Observatory. Narrow widths of synchrotron X-ray filaments [51, 52] and then- 
short time variability [53] require strong magnetic fields. Acceleration of protons up to the knee 
in the CR spectrum at ~ 10 15 eV can be realized only if magnetic field amplification occurs. 
Gamma-ray observations can put an important constraint on the strength of magnetic field, since 
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Figure 7: GeV luminosity (0.1-100 GeV) in units of 10 35 erg s“' plotted as a function of diameter squared for the SNRs 
detected with the Fermi LAT. Filled circles correspond to the SNRs interacting with molecular clouds, Rectangles are 
young SNRs. Pup A, middle-aged without clear molecular interactions, is shown as an open circle. 

the y-ray flux due to IC scattering is scaled simply by magnetic field strength given the observed 
flux of synchrotron radio emission; 

(iii) Acceleration spectrum— In the test-particle regime, the DSA theory predicts the spec- 
trum of CRs accelerated at SNR shocks obeying a power law Q(E) x E~ s with s x 2. The 
CR spectra measured directly at Earth have the form N{E) x E ~ 21 below the knee energy, as 
confirmed by the Fermi LAT measurements of the Galactic diffuse y-ray emission. Then the 
diffusion coefficient for CR propagation should be roughly D(E) x Q(E)/N(E) x £° 7 , Such a 
strong energy dependence would be inconsistent with the CR anisotropy measurements [54, 55]. 
Recently it has been suggested that the finite velocity of scattering centers (usually neglected 
in the DSA theory) due to magnetic field amplification may lead to steeper spectra of acceler- 
ated CRs with s ~ 2.2 [56], though the prediction depends on poorly understood properties of 
the magnetic turbulences produced at shocks. In this regard, measurements of the GeV y-ray 
spectrum, especially that of hadronic origin, are of great importance. 

62.1. Cos A 

Cas A is thought to be the remnant of a Type lib supernova that exploded around A.D. 1680. 
Since it is at the early Sedov-Taylor phase, the accelerated CR spectrum can reach already its 
maximum attainable energy and the process of energy conversion, from the expansion kinetic 
energy into CR energy, can well proceed. Cas A was the first SNR detected in the TeV y-ray 
band [57]. Though the emission mechanisms responsible for the TeV y rays remained unsettled, 
their detections provided direct evidence for particle acceleration beyond 10 TeV. Observations 
with EGRET resulted in setting upper limits on the y-ray flux above 100 MeV. 
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Figure 8: Gamma-ray spectra of young SNRs: Cas A (red) and RX 11713.7-3946 (black). Cas A: the I .AT spec- 
trum in 0.5-50 GeV (filled circles) (47J is shown along with the TeV spectra reported by HEGRA (line) [57], MAGIC 
(squares) [60], and VERITAS (open circles) [61]. RX J1713.7-3946: the LAT spectrum in 0.5-400 GeV (filled cir- 
cles) [48] is shown in combination with the energy spectrum of HJE.S.S. -detected emission (small filled circles) [62]. 
The y-ray spectrum of Cas A is best described by a hadronic emission model (n°-decay y rays), while the y-ray emission 
from RX J1713.7-3946 can be fitted by a leptonic model (inverse-Compton scattering). 


A clear detection of Cas A at a significance level of 12cr was made with the Fermi LAT 
using the y-ray data during the period between 2008 August 4 and 2009 September 4 [47], The 
GeV y-ray data were consistent with a point-like source, as expected from the angular size of 
the remnant (a radius of 2.5'), and also consistent with steady emission. The y-ray spectrum 
measured for Cas A in the GeV-TeV band is shown in Figure 8. 

The X-ray observations of synchrotron-emitting filaments measure the magnetic field at the 
forward shock; B ^ 0.3 raG is inferred from the observed width of the filaments [58], which 
is consistent with the variability timescale of the synchrotron filaments [59]. By combining the 
estimated magnetic field with the synchrotron spectrum in the radio bands, it was shown that the 
leptonic y-ray emission model does not well explain the LAT-detected emission 147]. Instead, the 
GeV-TeV y-ray emission can be ascribed mainly to the a 0 -decay component. The total proton 
content at the current age of the remnant amounts to W p » 0.4 x 10 50 erg for the shocked gas 
density of = 10 cm -3 , which is estimated from the measured dynamics of the remnant. This 
would be the first case where one can estimate the amount of CRs in SNRs with a reasonable 
accuracy. The spectral slope of accelerated CR protons is assumed to be s = 2.3 from the radio 
spectral index (Figure 8), which seems compatible with the GeV-TeV spectrum. Both the CR 
content and spectrum are in reasonable agreement with the direct measurements of the CR flux 
and the CR anisotropy. 

6.22. Tycho 

Tycho’s SNR is the remnant of a Type la explosion in 1572, currently in transition from the 
ejacta-dominated phase to the Sedov-Taylor phase. Recent detection of TeV y-rays from Tycho’s 
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SNR by the VERITAS Collaboration [63] has just been followed by the very recent 5 <r detection 
of GeV y-ray emission with the Fermi LAT [64]. The LAT spectrum is characterized by photon 
index T = 2.3 ± 0.2(stat)±0.1(sys), which smoothly connects with the VERITAS spectrum. 

Given the amplified magnetic field of 0.2-0. 3 mG estimated for shock downstream [65, 58] 
and the constraint on the ambient density (nu < 0.3 cm -3 ), the y- ray spectrum of Tycho’s SNR 
can be explained only by the 7r°- decay emission. Compared with Cas A, the ambient density is 
smaller and the radio emission is much weaker, so that the leptonic scenario faces a more serious 
difficulty. The total CR proton content at the current age amounts to W p * (0.6-1. 5) x 10 50 erg 
for a reasonable range of SNR parameters such as distance. The accelerated CR protons and 
electrons with number index s = 2.3 are inferred from the multiwavelength fitting. Such a steep 
spectrum may be explained by the finite velocity of the scattering waves in the shock upstream 
region [66]. As for Cas A, the LAT results obtained for Tycho’s SNR strengthen the case for 
SNR origin of the Galactic CRs. 

6.23. RXJ 17 13. 7 -3946 

SNR RX J17I3.7-3946 (G347 .3-0.5 in the Green’s SNR list [67]) is well known for its 
bright synchrotron X-ray emission that completely dominates the radiation output from the rem- 
nant [68, 69] and also for its bright TeV y-ray emission that well traces the synchrotron X-ray 
map [70, 62], Since its discovery, the origin of the TeV y-ray emission has attracted strong 
interest from the high-energy astrophysics community [71 , 72]. 

Figure 9 presents the Test Statistic (TS) map of a sky region around SNR RX J1 7 13. 7-3946 
calculated using the events above 0.5 GeV taken from 2008 August 4 to 2010 August 4 [48]. 
(TS is defined as TS = -21n(Lo/Lpomt) where Lq is the likelihood value of the null hypothesis 
that takes account of only background sources in the model, and Lp^, t is that of the hypothesis 
assuming a point source at each location.) The extended y-ray emission coincident with SNR 
RX J1713.7-3946 can be seen in Figure 9, though it was hard to investigate internal spatial 
structures due to the lack of photons. 

In Figure 8, the energy spectrum of SNR RX J1 713. 7-3946 measured with the Fermi LAT 

[48] is shown together with that of the H.E.S.S. telescopes [62], The LAT spectrum can be 
characterized by a hard power law with photon index T = 1.5 ± 0.1(stat)±0.1(sys), smoothly 
connecting with the steeper TeV spectrum. The hard power-law shape agrees with the expected 
IC spectrum (the leptonic model), as illustrated in Figure 8. Given the energy flux ratio of the 
observed synchrotron X-ray emission and the y-ray emission, the leptonic model requires that 
the average magnetic field be weak, B =* 10 /zG [62]. If the leptonic scenario is the case, the 
filamentary structures and variability in X-rays [53] would be attributable to locally enhanced 
magnetic fields, e.g., time-variable local turbulent magnetic fields [73], Importantly, within the 
leptonic scenario the y-ray spectrum provides a robust estimate of the total amount of relativistic 
electrons, as W e ~ 1 x 10 48 erg. 

6.3. SNRs Interacting with Molecular Clouds 

The brightest GeV y-ray sources associated with SNRs are middle-aged remnants that are 
interacting with molecular clouds (MCs) as evidenced by the detections of 1720 MHz OH maser 
emission [75], Figure 10 presents the raw LAT count maps in 2-10 GeV of the four SNRs (W51C 
[41], W44 [42], IC 443 [43], and W28 [44]) obtained with 2.5 yr data of the LAT observations 

[49] . Even before the subtraction of the intense Galactic diffuse emission, the y-ray signals from 
the SNRs can be easily seen. 
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Figure 9: Test Statistic maps for a point source hypothesis at each location of the sky in the region around SNR 
RX J1713.7-3946 [48], Events above 0.5 GeV are analyzed. TS maps were obtained by maximum likelihood fit- 
ting, where (a) the Galactic and isotropic diffuse emissions and 1FGL catalog sources [74] were included as background 
sources, (b) Sources A,B,and C were also taken into account in the background model. 



Figure 10: Fermi LAT 2-10 GeV count maps (in celestial coordinates) of the MC-interacting SNRs with extended y-ray 
emission: (a) W51C; (b) W44; (c) IC443; and (d) W28 [49]. The intensity is in units of counts per pixel with a pixel 
size of 0.05 X 0.05 deg 2 , Superposed are the contours from the VLA radio maps (see [49] and references therein). A 
black ellipse in panel (a) represents the location of shocked CO clumps. The black crosses in panel (c) are the locations 
of shocked molecular clumps from which OH maser emission is detected. The position of a pulsar wind nebula is also 
marked in panel (c). 
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Figure II: Gamma-ray spectra of shell-type SNRs measured with the Fermi LAT: W51C (filled circles [41]); W44 (open 
circles [42]); IC 443 (filled rectangles [43]); W28 (open rectangles [44]); Cassiopeia A (filled diamonds [47]). 

The MC-SNRs (namely, SNRs interacting with MCs) constitutes the dominant class of LAT- 
detected SNRs (Figure 7). The y-ray luminosity in the 1-10 GeV band spans Ly = (0.8-9) x 
10 35 erg s' 1 , larger than for the young SNRs like Cas A and RX J1713.7-3946 (Figure 1 1). The 
GeV y-ray spectra commonly exhibit a spectra] steepening in the Fermi LAT band, demonstrating 
the importance of LAT observations of these SNRs. Also, the synchrotron radio emission of the 
four SNRs can be characterized by a large flux density of 160-310Jy at 1 GHz with flat spectral 
index of a ^ 0.3-0.4 [67]. The LAT-detected SNRs are generally radio-bright objects [49]. 

The predominance of the MC-SNRs among the SNRs detected by the Fermi LAT and their 
high y-ray luminosity indicate that the y-ray emission should be enhanced by the interactions 
with molecular clouds. Nonthermal bremsstrahlung by relativistic electrons or 7r°-decay y- rays 
produced by high-energy protons are the two plausible channels of the y radiation. As discussed 
above, the observed high luminosity of the GeV y-ray emission disfavors the IC scattering off 
the CMB and interstellar radiation fields as the emission mechanism. 

There are two different types of models to account for the GeV y-ray emission from the 
MC-SNR systems. The Runaway CR model considers y-ray emission from molecular clouds il- 
luminated by runaway CRs that have escaped from their accelerators, namely SNRs [76, 77,78]. 
If the y-ray emission detected by the Fermi LAT can be ascribed to the escaping CRs, one can 
leam about how CRs are released into interstellar space and how they propagate in the vicinity 
of the SNR where self-generated Alfven waves may change the local diffusion coefficient [79]. 
Another scenario is the so-called Crushed Cloud model [80, 81] that invokes a “shocked” molec- 
ular cloud into which a radiative shock (typically v s ~ 100 km s' 1 ) is driven by the impact of 
the SNR’s blastwave. The cosmic-ray particles accelerated at a cloud shock are adiabatically 
compressed behind the shock front, resulting in enhanced synchrotron and 7r°-decay y-ray emis- 
sions. As shown in Figure 12, reacceleration of pre-existing CRs at a cloud shock is capable 
of explaining both radio synchrotron and GeV y-ray emission in the case of SNR W44 [81], 
The y-ray luminosity of Ly ~ 10 35 erg s' 1 in 1-10 GeV agrees with the theoretical expectation. 
Morphologies of synchrotron radio emission support this scenario. The radio continuum map of 
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Figure 12: Radio (left) and y-ray (right) spectra of SNR W44 together with the reacceleration model [81], 

SNR W44 exhibits filamentary and sheet-like structures of synchrotron radiation well correlated 
with the shocked H 2 emission [82].' One can learn about shock-acceleration processes in dense 
environments, which complements the studies of young SNRs. 

Both sites of the -y-ray emission could be important. For example, a hard GeV y-ray source 
[44] which is also bright in the TeV band [83] outside the southern boundary of SNR W28 might 
be explained by runaway CRs, while 7r°-decay emission from radiatively-compressed layers of 
clouds could account for the GeV y rays observed within SNR W28. Detailed analysis of spatial 
distributions of GeV y rays combined with CO observations will make it possible to discern the 
y-ray sources that are attributable to runaway CRs. 

6.4. Other Candidate Cosmic-Ray Sources 

Although SNR remain primary candidates for accelerating the bulk of cosmic rays, other 
Galactic sources observed by the Fermi LAT also represent sites of particle acceleration and 
interaction. Such sources include pulsar wind nebulae [84], high-mass binary systems such as 
LS I +61°303 [85] and Cygnus X-3 [86], whose spectra seem most consistent with leptonic 
sources , and Nova V407 Cygni [87] , a symbiotic binary system whose spectrum can be modeled 
with either leptonic or hadronic processes. While the energetics requirement for the CR protons 
is generally difficult to satisfy by these non-SNR sources, they might be the dominant sources 
of the CR electrons and positrons. The large number of unidentified sources found in the Fermi 
LAT 1FGL catalog [74] offers the possibility that there may be yet-undiscovered CR production 
sources. 

7. Conclusions 

Both electron and y-ray observations made with the Fermi Large Area Telescope offer new 
insights into cosmic-ray acceleration and interactions. The direct electron plus positron measure- 
ments suggest, although they do not prove, a possible nearby source. Future LAT observations 
will help measure separately the positive and negative components up to 200 GeV and will im- 
prove on the search for anisotropies already begun. On a Galactic scale, y rays suggest a larger 
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scale height or flatter radial distribution of cosmic rays than previously expected. From the y- 
ray studies of supernova remnants, LAT has found growing evidence that at least some SNR are 
likely accelerators of hadronic cosmic rays, although one example, RX J1713.7— 3946, appears 
more consistent with being a leptonic source. The Fermi LAT continues to collect data. Improved 
analysis methods, coupled with the higher statistics from the ongoing mission, can confidently be 
expected to produce even more results addressing questions about cosmic-ray origin and propa- 
gation in the future. 
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